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ABSTRACT 

The Sun's corona is millions of degrees hotter than its 5,000 K photosphere. This 
heating enigma is typically addressed by invoking the deposition at coronal heights of 
non-thermal energy generated by the interplay between convection and magnetic field 
near the photosphere. However, it remains unclear how and where coronal heating 
occurs and how the corona is filled with hot plasma. We show that energy deposition at 
coronal heights cannot be the only source of coronal heating, by revealing a significant 
coronal mass supply mechanism that is driven from below, in the chromosphere. We 
quantify the asymmetry of spectral lines observed with Hinode and SOHO and identify 
faint but ubiquitous upflows with velocities that are similar (50-100 km/s) across a 
wide range of magnetic field configurations and for temperatures from 100,000 to 
several million degrees. These upflows are spatio-temporally correlated with and have 
similar upward velocities as recently discovered, cool (10,000 K) chromospheric jets 
or (type II) spicules. We find these upflows to be pervasive and universal. Order 
of magnitude estimates constrained by conservation of mass and observed emission 
measures indicate that the mass supplied by these spicules can play a significant role in 
supplying the corona with hot plasma. The properties of these events are incompatible 
with coronal loop models that only include nanoflares at coronal heights. Our results 
suggest that a significant part of the heating and energizing of the corona occurs at 
chromospheric heights, in association with chromospheric jets. 

Subject headings: Sun: atmospheric motions — Sun: magnetic fields — Sun: chro- 
mosphere — Sun: transition region — Sun: corona 
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1. Introduction 



Statistical correlations between the brightness of the chromosphere and corona strongly sug- 
gest that these two regions are powered by a similar driver dSchrijver & ZwaanlllOOOl) . Yet, despite 
some suggestions that much of the non-thermal energy driving coronal heating is depos ited in 
the chromosphere (|Aschwanden et alj 120061 : iGudiksen & Nordlundl 120051: lMcIntoshll2007h . most 
models of coronal loops treat the chromosphere as a passive mass reservoir that merely reacts to 
energy release at coronal heights and t he associated cha nges in the coronal thermal conductive flux 
or deposition of non-thermal energy (|Klimchukl |2006|) . The differing role of the chromosphere 
in these models implies that observations of how the corona is filled with plasma from the chro- 
mospher e hold significant diagnostic potential for d istinguishing between various coronal heating 
models ( Klimchukll2006l : |Patsourakos & KlimchukI 12006) . Yet, it has been very difficult to estab- 
lish how the dynamics and energetics of chrornosphere and corona are cou pled (iKlimchukI 12006 : 
Mcintosh et allboOTi : kiouaz. Peter & Lemair^boOTl : iHansteen et allboovk 



Chromospheric spicules have long been considered a candidate for such coup l ing by providing 
mass to the corona and/or solar wind feeckers 1968 : Pneuman & Kopp 1978 : Athay & Holzer 
19821). These jet-like features, with lifetimes of order 3-10 m inutes, propel cool matter upwards 
to coronal heights with speeds of 20-30 km/s (|Beckerslll968|) . A signifi cant fraction of these are 
likely caused by magnetoacoustic shocks on rn agnetic flux concentrations ( De Pontieu et al.l2004a : 
Hansteen et al.ll2006l : lOe Pontieu et al.ll2007a|) . Spicules are observed at temperatures from 5,000 
to 500 ^,000 K and are estimated to carry upward a mass flux 100 times larger than that of the solar 
wind (jPneuman & Kopplll978|) . However, no s i gnatures of the se ejecta at coronal temperatures 
have been reported previously (|Withbroelll983l : iMariskal 119921) . Consequ ently, a direct r ole of 
spicules in the coronal mass/energy balance has been dismissed as unlikely (|Withbroelll983|) . 



However, the high spatio-temporal res olution of the Solar Optical Telescope (SOT: lTsuneta et al. 



2008|) on board Hinode (iKosugi et al.ll2007b has revealed a previously unrecognized second type of 
spicules (iDe Pontieu et al.ll2007bl) that is a prime candidate for establishing a link between corona 
and chromosphere. These "type 11" spicules have much shorter lifetimes (10-100 s), are more 
violent (upflows of order 50-150 km/s), and often fade rapidly from the chromospheric Ca II H 
3968A SOT passband, which suggests rapid heating and ionization of singly to at least doubly ion- 
ized calcium. The formation of these ubiquitous type II spicules has been attributed to magnetic 
reconnection ( De Pontieu et al.ll2007bl) . 



To study the thermal evolution of these spicules, we use simultaneous observations of Ca II 
H emission (with SOT) and EUV spectral lines formed in the corona at temperatures of s everal 
million K from the Extre me-ultraviolet Imaging Spectrome ter (EIS: ICulhane et al.ll2007|) . We 
use spectra from SUMER dwilhelm et al.lll995l) on SOHO dFleck et al.lll995l) to show that the 
process we describe is ubiquitous and occurs in active regions, quiet Sun and coronal holes. The 
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observations and our novel analysis techniques are described in §[2l while §|3]covers the results of 
our analysis. In §|4]we discuss the impact of our results on the existing coronal heating hypotheses. 



2. Observations & Analysis 

We focus on plage regions in active regions NOAA 10976 and 10977 observed near disk 
center on 2 and 5 December 2007 by Hinode/SOT and EIS. By observing at disk center we reduce 
the ambiguities introduced by line-of-sight superposition at th e limb (where sp icules are most 



readily identified), which has plagued such studies in the past (|Withbrodll983[) . SOT obtained 
a timeseries of images using the Ca II H 3968 A filter (with 0.5s exposure and 8s cadence) that 
covered a 25"x40" area of plage. Simultaneous EIS observations of the same region started with 
a region-wide scan with coarse horizontal (5") stepping with 60s exposures for a wide variety of 
spectral lines. This was followed by much finer (1" stepping), rapid (60s exposure) coronal rasters 
with EIS over a 5"x80" sliver of plage. 

To connect the chromospheric and coronal dynamics we need to coalign the coronal spec- 
tral rasters with the high cadence chromospheric imaging. This is not straightforward. The SOT 
timeseries are coaligned onboard Hinode using a correlation tracker. The EIS spectra are pointed 
independently from SOT, and undergo significant jitter caused by the spacecraft, and thermal flex- 
ing within EIS. We remove the spacecraft jitter from the EIS pointing by using the xrtjitter IDL 
routine (from solarsoft). To correct for the effect of thermal flexing on EIS pointing, we use sev- 
eral steps. First, we perform additional coalignment of the Ca II timeseries using cross-correlation. 
This removes drifts caused by solar evolution within the small field of view of the SOT correlation 



tracker. We use chromospheric and coronal images from TRACE (|Handy et al.lll999|) to connect 



the EIS and SOT data and determine the EIS pointing as a function of time. 

This pointing information is used to calculate synthetic rasters of upper chromospheric ac- 
tivity (UCA) derived from the Ca IIH timeseries. The latter is dominated by slowly evolving (on 
timescales > 60s) photospheric and low-chromospheric contributions because of the wide band- 
pass (FWHM=2.2A) of the SOT Ca II H filter. To isolate the relatively faint signature (from the 
Ca H core) of the upper-chromospheric type-II spicules in the Ca II H ti meseries, we exploit their 



highly dynamic nature: they evolve on much shorter timescales (< 60s) (iDe Pontieu et al.ll2007bl) 
than the photospheric and low-chromospheric contributions. We thus perform, for each pixel, tem- 
poral high-pass Fourier filtering on the timeseries with a cutoff of 18 mHz. The resulting timeseries 
of (the absolute value of the) high frequency signal is a good proxy for UCA because it removes 
the photospheric and low-chromospheric features (online movie 1). 



To study the high temperature response, we analyze EIS spectra of moss regions which 
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constitute the transition region (TR) at the footpoints of active region loops ( 


Berger et al. 


1999; 


De Pontieu et al. 


1999: 


Fletcher & De Pontieu 


1999; 


Martens. Kankelborg & Bergei 


2000 


). We 



find that the bright core of the coronal spectral lines are typically Doppler shifted by less than 
10-20 km/s and note that the Doppler shift and intensities of th e coronal profiles sho w no spatial 
correlation betwe en upper-chromospheric and coronal activity (IHansteen et al.ll2007h . However, 
Hara et al.l (|2008|) recently discovered that there is often a faint blue-shifted component in addition 
to the bright core of the line, indicating coronal up flows of order 50-100 km/s. To isolate these 
upflows from the bright core of the spectral line more clearly, we calculate maps of the blue-red 
(B-R) asymmetry in the coronal spectral lines (Fig.fl]) and find that the upflows occur preferentially 
in moss footpoint regions. 

To calculate the B-R asymmetry of spectral lines, we first fit a single Gaussian to the line 
profile Ix, and determine the line centroid Aq of the Gaussian fit. The B-R asymmetry for an offset 
AAi from the line centroid is then given by: 



BR 



AAi 



Ao-AAi+5A„ 
Ao — AAi —5X^1 



Ao+AAi-5A^ 



(1) 



with SXu, the wavelength range over which the B-R asymmetry is determined. To determine the 
propagated errors on the B-R asymmetry measure, we use estimates for the error on the intensity 
using the EIS software (including photon noise and other errors). In our maps we only show 
locations where the B-R asymmetry is larger than twice the estimated propagated error. We apply 
the same method to SUMER rasters of qui et Sun and coronal h ole on 6-nov-1999, after calibrating 
the data using the techniques described by lDavey et al.l (|2006|) . 



3. Results 

We calculate B-R asymmetry maps for EIS rasters of Fe XIV 274A (formed at 2 MK under 
equilibrium conditions) for velocities arond 80 km/s and 120 km/s (Fig.|2j5,c). Both B-R asymme- 
try maps show no obvious correlation with the intensity of the core of the Fe XIV line (Fig.|2^) from 
which they are derived. In contrast, we find that the sites of faint upflows at coronal temperatures 
(Fig. |2j5,c) are well correlated with upper-chromospheric activity (Fig. [2jl, online movie 2). The 
correspondence between the B-R and UCA maps is surprisingly good, given several complicating 
factors. 

First, despite our careful coalignment, there remains uncertainty of order 2" in the EIS/SOT 
coalignment. In addition, the correlation between the chromospheric signal and coronal blueshifts 
associated with spicules depends greatly on the viewing angle between the line-of-sight and the 
direction of the upward flows (Fig. 5). When looking down at loop footpoints, we observe the 
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chromo spheric jet as a Ca II H brightening (because of projection), associated with a faint, strongly 
blue-shifted component in the coronal line. When the line-of-sight is not aligned with the upflows, 
the spicules will appear as jets in the chromo spheric line, but we no longer see blue-shifts in the 
coronal line since the line-of-sight velocity is reduced, and the spicular emission measure is too 
small compared to the dominant coronal emission. The wide range of angles between upflows 
(i.e., magnetic field) and line-of-sight thus naturally leads to an imperfect correlation between 
chromospheric brightness and coronal upflows. It is also possible that some of the UCA sites that 
do not match coronal upflow events do have coronal counter parts, but at much higher or lower 
temperatures than the Fe XIV line we observe. We also cannot expect a one-to-one correlation 
between the upper chromospheric brightness and the brightness in the coronal asymmetry maps, 
because we do not know the exact formation and heating mechanism of the spicules. Finally, some 
of the coronal upflow events could be nanoflares that are driven by energy deposition in the corona, 
i.e., without significant chromospheric brightness enhancement. 

The relationship between these chromospheric events and coronal upflows is elucidated by 
our finding that the velocities of the latter in two Fe XIV lines (Fig. [3]) are similar to the upward 



velocities for type II spicules seen at the limb (from lDe Pontieu et al.ll2007bl) . both peaking between 
50 and lOOkm/s. We find that these upflows are pervasive and universal: UV lines observed with 
SUMER show faint upflowing components with velocities of the same magnitude (50-lOOkm/s) 
that are well correlated across the TR and low corona (100,000 to 600,000 K) for a whole range 
of magnetic field configurations (active region plage, and network in quiet Sun and coronal holes, 
with the velocities in the latter two somewhat lower than in plage). 

This is illustrated with spectroheliograms of C IV 1548 A and Ne VIII 770A from for a coronal 
hole region near disk center (Fig. H]). We find a clear difference in Ne VIII intensity (reduced) and 
Doppler velocity (blueward) as well as C IV Doppler velocity (blueward) in the coronal hole 
compared to the quiet Sun. These differences can be understood in terms of a picture where the 
bulk of the material re aching the formation temperature of Ne VIII is injected into the solar wind 



(|McIntosh et al.ll2007|) . We also find that the average profile in the network in both quiet Sun and 
coronal hole shows significant asymmetry towards the blue for velocities of 50-100 km/s for both 
C IV and Ne VIII lines. The B-R asymmetry maps show that faint upflows at 50-100 km/s are 
well correlated for a lower TR (C IV) and a coronal (Ne VIII) line and preferentially occur in and 
around the network. The spatial patterning and correlation of these upflows does not change as one 
crosses the boundary between quiet Sun and coronal hole. This is surprising, since the magnetic 
field configuration in the corona is dramatically different. 
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4. Discussion 

Our analysis of the asymmetry of spectral lines is very sensitive to the presence of blends in 
the wings of the lines, which can skew the asymmetry maps towards the red or blue so that these 
are no longer indicative of flows. For that reason, we have avoided lines that have strong blends. 
We have extensively studied the influence of blends on the lines used, and find that none of them 
change our conclusions (see a follow-up paper for a detailed discussion). We point out that the 
upflow maps for different Fe XIV lines (at 264 and 274A) are very similar, and that we observe 
similar upflows for a wide range of temperatures (C IV 1548 A, Ne VIII 770A, Fe XIV 264 and 
274A). 

We estimate that these spicule-associated upflows play a significant role in replenishing the 
corona with hot mass. To play such a role in active regions, these hot upflows must provide enough 
mass flux = piViNi to replenish the corona which loses mass by downflows (after cooling) 
at a rate of /c = 'jPchc/tc- Here, pi is the density of the coronal mass propelled upward and 
heated in a spicule, Vi the upward velocity of the spicule, A^j the number of spicules that occur 
at any location over the spicule lifetime U, pc the density in coronal loops, he the coronal scale 
height, tc the time for coronal plasma to cool to chromospheric temperatures, and 7 the ratio of 
cross sections of the upper corona and the plage region below. Our observations in coronal lines 
also show that the faint blue-shifted spicular emission measure = {eTPiYvitiNi (with the 
fraction of total spicular coronal density at a temperature T) can only be a fraction a of that of 
the dominant emission Cc = p^hc at the core of the lines, which is emitted by the slowly cooling 
plasma in previously filled loops. We thus require that fi = fc and Cj = aCc, and find that 
Ni = '^'^e^h rtj/iavitl), and ft = aPrt rJ iie^tj). For realistic estimates of 7 = 3, /ic = 50, 000km, 



= 2000s (ISchriiver & Zwaanll2000h. a = 0.05 (from analysis of EIS spectra. Fig. [I]:,d), U 



100s and Vi = lOOkm/s (IDe Pontieu et al.ll2007bh we find ^ ^ Pc/(3e|,) and A^^; ^ 2e|,. Fe XIV 



observations only cover temperatures around 2 MK, and most likely reveal only a fraction of the 
coronal component of spicules, which may range in temperature from 1 to lOMK. If we make the 
reasonable assumption that only a fraction = 0.25 — 0.5 of the spicular coronal mass flux is 
observed in Fe XIV, we find, for coronal densities pc ~ lO^cm^'^, a total density of the coronal 
component of spicules of order 1 — 5 x lO^cm^'^. 

This implies that these spicules can fill the corona with hot plasma even i f only between 



1 and 5 % of the chro mospheric spicule densities reach coronal temperatures (IBeckersI Il968 
Pneuman & Kopplll978|) . We can also estimate how many spicules Np are required at any time 



in a plage region of diameter 6: Np = Ni5 /6f with 6i the diameter of a spicule. For 5 ~ 8000km 
and 5i ~ 250km, Np is of order 80-320, implying a spicule number density of 0.7-3 spicules per 
square arcsecond of plage. This density is compatible with SOT observations, which implies that 
these spicules can play a significant role in filling coronal loops. Individual spicular events have 
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energies of the order of nanoflares (~ 10^^ erg) and must recur every ti/Ni seconds at the same 
location (i.e, 250 -1000 s for the above estimates) to replenish the coronal mass. The energy flux 
carried into the corona by these hot spicules is of order 5 x 10^ erg cm~^ s^^ with the above as- 
sumptions and if we assume that all spicules reach 2 MK temperatures. This is of the same order 
of magnitude as the energy flux required to heat the active region corona (|Priestill982|) . 



The spatio-temporal correlation between upper chromospheric activity (associated with jets) 
and faint upflows at several MK, and the ubiquity of faint and correlated upflows with velocities 
of order 50-100 km/s for temperatures ranging from 10,000 to 2 MK for a variety of magnetic 
field configurations (active region, quiet Sun and coronal hole) poses significant challenges to our 
current understanding of the mass cycling of the corona. The properties of these events do not 
seem compatible with those predicted by coronal loop models that are driven only by nanoflares 
at coronal heights. Such models predict high-speed upflows at very high temperatures of several 
million deg rees, with either lower velocitie s or no detectable upflows at lower temperatures (<1 
million K) (|Patsourakos & Klimchukll2006|) . More generally, our observations of similar up flow 
velocities for a wide range of temperatures from 10,000 K to 2 million K are a challenge to any 
model in which upflows are driven by transient overpressure at chromosphe ric heights (e.g., by 
deposition of non-thermal energy from electron beams generated in the corona: iRaftery et al. I I2OO9I) 
because those lead to increasing velocities with temperature. 

Our finding that these upflows do not change qualitatively or quantitatively between quiet 
Sun and coronal holes strongly suggests that they are driven from below. This is because models 
in which energy is deposited at coronal heights are dependent on conductive flux and/or reconnec- 
tion at tangential discontinuities, which are expected to have significantly different properties in 
open versus closed field environments (|Parkeiill989|) . Such models also do not provide a natural 
explanation for why the chromospheric counterparts show upward extrusions and can be seen in 
lines from neutral elements such as Mg I b 5172 A (SOT). 

Nanoflares at coronal heights may be important for coronal heating, but they cannot explain 
our observations. Our analysis suggests a scenario (Fig. [5]) where a significant, potentially domi- 
nant, part of the heating of the coronal plasma occurs at c hromospheric heights, in association with 
chrom ospheric jets (similar to previous suggestions, e.g. JPneuman & Kopplll978l : lAthay & Holzer 
1982|) . A significant fraction of the plasma propelled upward in these chromospheric jets or type 
II spicules is heated to coronal temperatures, providing the corona with hot plasma. The coro- 
nal emission of individual spicules is faint compared to the dominant emission, which originates 
in previously filled loops that are slowly cooling. Previous efforts to establish a connection be- 
tween chromosphere and corona have focused on properties of the core of coronal lines. The 
weak emission of the coronal upflows associated with spicules, which prov ide the direct link to th e 
chromosphere, helps explain why these efforts have failed in the past (e.g. jHansteen et al.ll2007() . 
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Conceptually, our findings are compatib le with results from numerical simulations (IGudiksen & Nordlund 



20051 : iHansteen. Carlsson & Gudiksenll2007|) that suggest that most of the heating in the solar at- 



mosphere occurs below heights of 2 Mm, i.e., at chromospheric heights. Preliminary results of 
more recent simulations suggest a scenario in which magnetic reconnection in the chromosphere 
causes upward jets of plasma at a range of temperatures owing to magnetic Lorentz forces even as 
the plasma is heated to TR and coronal temperatures in the process. 
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Fig. 1.— Fe IX/X 171 A images (from TRACE) and EIS rasters in Fe XIV 21 Ak of NOAA AR 
10977 on 5-dec-2007 show cooler moss footpoints (log T=5.9, a) as well as hotter loops (log T= 
6.2, b) connecting the footpoints. Spectral line profiles of the coronal Fe XIV line show a deviation 
from a Gaussian (dashed line, c) indicative of hot plasma flowing upwards at high speeds. Spatial 
maps of the B-R asymmetry at 50-120 km/s (d, see text for details) show a good correspondence of 
blueshifted events with the locations of moss footpoint regions (a). Online movie I of the region 
indicated by the white box in panel (a) shows how high-pass Fourier filtering isolates the upper- 
chromospheric activity (right panel of movie) from the more slowly evolving photospheric signals 
that dominate the Hinode Ca II H bandpass (left panel of movie). 



- 12- 




Fig. 2. — The locations of blueshifted signals from the B-R asymmetry maps of Fe XIV 274A 
(panel a shows the intensity) centered around 120 (panel b) and 80km/s (panel c) correlate well 
with a synthetic raster map of the upper chromospheric activity (panel d, which is derived from 
panel e) of a small plage region on 2-dec-2007. Online movie 2 shows a timeseries of similar 
comparisons. 
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Fig. 3. — A histogram of upward velocities in chromo spheric type n spicules (black line) as 



derived from Ca II H 3968A timeseries at the limb (from lDe Pontieu et al.ll2007b|) shows a range 



of velocities that is similar to coronal upflow velocities derived from subtracting a Gaussian fit 
from the profile of Fe XIV 264A (blue) and Fe XIV 274A (green) lines for all locations that show 
significant upflows in active region 10977 on 5-dec-2007. 
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Fig. 4. — Panels a and b show the logarithm of intensities of C IV 1548 A and Ne VIII 770A for 
an equatorial hole observed at disk center with SUMER. The solid contours are a proxy for the 
supergranular network boundaries (derived from a threshold of 150 counts in panel a). The solid 
white contour delineates the boundary of the coronal hole from the quiet sun. Panels c and d show 
the Doppler shift of the lines. Below this we show representative normalized network line profiles 
taken from x=-300",y=-42". We determine B-R asymmetry maps (panels e and f) for a 24 km/s 
wide region [shown as red/blue rectangles] at offsets from line center of 50km/s for C IV and 
60km/s for Ne VIII. In both cases, we see strong enhancement of the blue wing emission in the 
magnetic network. 
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Fig. 5. — Cartoon illustrating the mass and energy transport between the chromosphere, TR and 
corona, as deduced from SOT and EIS observations. See §4 for details. 



